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Abstract. Using the Mathematica package, we find exact analytical expressions for the 
so-called de-projected de Vaucouleurs and Sersic laws as well as for related spatial (3D) 
quantities - such the mass, gravitational potential, the total energy and the central 
velocity dispersion - generally involved in astronomical calculations expressed in terms 
of the Meijer G functions. 
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1. Introduction 

Dynamical studies of astronomical systems like Elliptical Galaxies or Clusters of Galaxies 
involve the de-projection of observed (projected on the sky) quantities like surface bright¬ 
ness profiles, numerical density profiles, velocity dispersion profiles etc.The 3D profiles 
obtained are then used to derive e.g the total luminosity (or mass) of the system or the 
gravitational potential and are used in the Jeans equation which is then resolved to get 
for instance the kinematics of the system. 

The de Vaucouleurs profile (de Vaucouleurs 1948) and its generalization by the Sersic 
law (Sersic 1968), is one of the most often used laws particularly in the study of Elliptical 
Galaxies. 

Unhappily these laws have so far lead to non-analytical de-projection (i.e spatial) 
quantities. Efforts have been made in the last decades to provide either numerical tables 
(Poveda et al 1960, Young 1976) or approximations and asymptotic expressions (Mel- 
lier and Mathez 1987, Ciotti 1991, Graham and Colless 1997, Ciotti and Berlin, 1999, 
Marquez et al. 2001). 

Here we give analytical exact expressions for 3D quantities usually derived when using 
the de Vaucouleurs or Sersic laws. 

2. Principles 

The classical de Vaucouleurs and Sersic laws express the dependence on the projected 
central distance R of for instance the Luminosity Intensity / (R) of an Elliptical Galaxy. 

The de Vaucouleurs law relates the Intensity I{R) to the central one, /(O), by: 
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= exp -7.66925 [[R/Re)^ - 1] 


( 1 ) 


and the Sersic profile, which is a generalization of the de Vaucouleurs profile, is 
written: 


I{R)/Io = exp -b{m) - 1] 


( 2 ) 


The parameter 6(m) is determined from the definition of the effective radius Re, which 
is the projected radius inside which the projected luminosity (or mass) equals half of the 
total luminosity (or mass). The de Vaucouleurs law is recovered for 6(4) = 7.66925. 

The requested 3D profiles are related to the derivative of the projected profiles by the 
usual Abel Integral written here for the 3D density profile n(r): 


i{r) = -- 


dl 


1 


,dR 


(3) 


TT Jr dR 

Except for some particular cases there is no known exact expression for theses inte¬ 
grals, in particular in the case of the Sersic ( de Vaucouleurs ) profiles. 

However, using Mathematica we succeeded in obtaining exact analytical expressions 
for such integrals that involve the Meijer G functions. Physical quantities such as the 
spatial luminosity or mass profiles, the gravitational force, the gravitational potential and 
energy, which are combinations or integrals of the above functions, have also analytical 
expressions involving Meijer G functions. Some of these expressions are given below and 
numerical evaluations are compared to previous numerical calculations. 

We first give the classical definition of the Meijer G functions together with some of 
their properties which will be useful in understanding the results given by Mathematica. 


3. The Meijer functions 


These functions are defined as integrals of products of T functions. The generalized Meijer 
G function is defined as: 

r f [r(l — oi — ru)... r(l — ai — ru)T{l3i + ru)... T{Pk + ru)] 


z ^du (4) 


2TTi J [r(Q:/+i -I- ru). .. r(a p-|- ru)r(l - Pk +i - ru)... r(l - Jq - ru)] 

(Gradshteyn & Ryzhik, 1980; Wolfram, 1991; for a collection of formulae related 

to the Meijer G functions cf. |ittp: //functions.wolfram.com/Hypergeometric Func- 

tions/MeijerG/). 


The case r = 1 defines the standard Meijer function, Gp’g 
which is the form we will be dealing with in the rest of this work. In the Mathematica 
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StandardForm notation it writes as: 




MeijerG {{ai ••• az}{Q:/+i • • • ap}}{{/3i • • •/3fc}{/3fc+i • • •/J,}},z = 

{oi • • • a/}, {ai+i ■ • • ap}\ 

{/3i • • • /3fc}, {fik+i ■ ■ ■ fiq} j 

For clarity we will keep both these notations, although suppressing the suffix “Meijer”. 
Notice that in these formulae, the empty case: { } means that the corresponding coeffi¬ 
cients are not defined and thus do not exist. This may occur, for instance, when k = q, 
01 1 = p, or else when one of those indices are null. An example is given by Eq. (|^) below. 

The following identity may be easily obtained by a substitution of the integration 
variable u ^ u + c, where c is a constant, in Eq. i): 






i-M) 




|{ai-|-C"ai-|-c},{ai4.i-|-C'"C(p-|-c}\ 


(5) 


{/3i-t-c---/3fc-t-c},{/3fc+i-t-c---/3g-t-c} j 

The moments of the Meijer function are expressible in terms of the higher-order Meijer 
functions: 




\{l3i---0k} ,{l3k+ 


i (jz — (2 


},{ai + i"ap} 


( 6 ) 


This may be straightforwardly demonstrated by inverting the order of the integrations. 

In some particular cases the Meijer G functions may be expressed in term of more 
classical special functions. As an example we give below the case of the Meijer function 


Gq ’2 (z|/3i,/ 32) (c.f. http://functions.wolfram, com/07.09.03.0330): 


G, 


(z|/3i,/32) = Gl% 

where Kr{x) is the modified Bessel function of r— order. 


(7) 


4. The 3D laws 


In the following, we give the analytical expressions using the Matheinatica formalism for 
the G Functions. 

Let first start by introducing some useful dimensionless quantities. The dimensionless 

2D and 3D x and s radial distances are expressed in terms of Re as: 

X = R/Re ; s = r/Re (8) 

The dimensionless 2D and 3D profiles are also defined as: 

i{x)=I{R)/Io and v{s) = n{r)^ (9) 

lo 

which lets equations (0) and transform in the following reduced form: 


i{x) = exp —b(m)(x^ — 1) 


( 10 ) 


where the de Vaucouleurs law is obtained for m = 4. 
The de-projection integral p) also transforms to: 


1 


v{s) = - 

IT 


di 


1 


dx Va;2 - 


dx 


( 11 ) 
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4.I. De-projection of the de Vaucouleurs law 


We give here some detailed results for the case of the de Vaucouleurs law as an example 
and give more general results in the next section using the Sersic law. Integrating the 
preceding equations formally using Mathematica, we first derive the expression for the 
3D profile n{r). We will then calculate the luminosity (or mass) profiles as well as the 
gravitational potential and the gravitational energy. 

Only numerical estimations or asymptotic behaviors were given before (Poveda et al, 
1960; Young, 1976; Mellier & Mathez, 1987) so we will compare our results with those 
provided by Young for the spatial density and for the luminosity (or mass). 

The basic ingredient to obtain the zz(s) profile is the derivative of the de Vaucouleurs 
law, equation (|^), relative to the projected dimensionless distance, x, which is written: 

=-^ exp [—b{x'^ — 1)] with b = 7.66925 (12) 

dx 4x4 

Integrating equation ( 0 ), the spatial density z^(s) expressed in terms of the dimen¬ 
sionless 3D radial distance s is then given by: 

zz(s) = 6.23828s-^/^ G({{}, {}}, {{i, || I 1,1, ^}, {}}, 0.713351 s^) (13) 


One then obtains the luminosity (or mass) spatial profile defined by: 

M(s) =47r f s''^v{s')ds' (14) 

Jo 

which with a new formal integration gives: 

Mis) = 39.1962 si G({{^}, {}}, {{i, 1,1,1, ^}, {- }}, 0.713351 s^) (15) 


The gravitational potential is defined by: 

which gives: 

«'(s) = 19.5981 si * 

II 111 5 3 7 7 ^ 9 5 

ltt8’8’4’8’8’ ’8’4-'’^ 


(16) 


}}, 0.713351 s^ 


(17) 


The values tabulated by Young (1976) are recovered defining the new potential: 

4''(s) = 1 — 4>(s)/4'(oo) 


such that 

4''(0) = 1 and 4>'(oo)) = 0 


In table 1, we compare the values given with the above expressions (normalized by 
the factor exp(6)(7r8!/6®)) to the numerical values obtained by Young for s ranging from 
10-® to 10. 
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Table 1. The spatial density and mass profiles for the de Vaucouleurs profile 


s 

u{s) 

i^Youngi^s') 

M(s) 

^Young (s) 

MO"® 

2.5553-10® 

2.5535-10® 

1.4730-10"“ 

1.4717-10"'' 

I-IO"® 

3.5797-10® 

3.5786-10® 

2.1130-10"''’ 

2.1122-10"® 

1-10"'‘ 

4.2189-10"' 

4.2183-10"' 

2.5959-10"'^ 

2.5954-10"’' 

I-IO"® 

3.7044-10® 

3.7042-10® 

2.4545-10"® 

2.4544-10"® 

MO"^ 

1.9679-10® 

1.9679-10® 

1.4961-10"® 

1.4960-10"® 

I-IO"^ 

4.4047 

4.4047 

4.4102-10"® 

4.4102-10"® 

1. 

2.1943-10"® 

2.1943-10"® 

4.1536-10"' 

4.1536-10"' 

10. 

7.8166-10"® 

7.8165-10"® 

9.4308-10"' 

9.4308-10"' 


J^.2. The de-projection of the Sersic law 


We give now more general expressions for the 3D profile derived from the Sersic law. 
This law is parametrized by an index m and a parameter b{m). The Sersic law for m = 1 
corresponds to a 3-D Exponential profile, already discussed by Fuchs & Materne (1982) 
and for m = 4 to the usual de Vaucouleurs law, treated in Section 2 above. In general. 


b{m) is found as a solution of equation (A.4) given in the Appendix. Values of b{m), (as 
well as logL(m)) have also been given before by Ciotti, Graham & Colless and Ciotti & 
Bertin. We give in the Appendix, values of b{m) calculated using Mathematica for the 
range m = 1 to m = 15 and compared to those of Ciotti & Bertin. 


The spatial density expressed in terms of the dimensionless 3D distance s can be 
written as (see e.g. Ciotti, 1991; Graham & Colless, 1997): 


1 


zz(s) = —/^'exp(-^)-p_ 

^ ^ TT Jo t ^ y/t-^rn _ ^ 

where b means 6(m) 


dt 


(18) 


Again this integral is found in terms of Meijer G Functions which reduces to a Bessel 
function (Kq) for to = 1. We give below the various expressions for various values of m 
and b{m). Defining auxiliary constants: 

, ( 19 ) 

(2TT)™‘y/m 

and 

( 20 ) 


( b{m) 
C2[m) = ' 


2 to 


we obtain, for to = 1: 
j^i(s) = 2ci{l)KQ{b{l) s) 


(21) 
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For m = 2: 

2.58697s-i G({{}, {}}, {{0, ^ ^ i}, {}}, 0.710231 


( 22 ) 


For m = 4: 

5.26885s-3G({{},{}},{{0 


113 3 15 3 
’8’4’8’8’2’8’4 


},{}}, 0.713351 


(23) 


Notice that Eq. (^^ is equivalent to Eq. dil) for the de Vaucouleurs profile given before. 
This may be seen by applying the identity Eq. (H) and noting that Ps = Pdv ~ 1/2, 
where Ps and Pdv denote the arrays of P coefficients appearing respectively in Eq. (^) 
and Eq. (|l^). 

The expressions above have the general form: 

Vm{.s) = Ci(TO)s^G({{},{}},{{/3s(ni)},{}},C2(TO)s^) 

= ci(m)s^Go™/? J) 

with Ps{m) denoting the 2m-array: 

Ps{m)^{ I 

( \ 2m J V 2m J m+l<j<2m J 


(24) 

(25) 


4-3. Other related quantities 

Related quantities as the mass , the gravitational potential, the (total) potential energy 
and the (central) velocity dispersion, can then be formally calculated by other integra¬ 
tions, similar to what has been done for the de Vaucouleurs profile. From Eq. ( p^ we 
find for the mass: 


M(s) = 27rci(m)s"™^'G ({{-^},{}},{{/3s(m)},{-^^^^}},C2(m)s^ 


^ 2^ci(m)s^G?:^;;:+i ( C2(m)s^ ^ ^ M (26) 

Notice that in the case m = 1, because of Eq. (il), an alternative expression for the 
mass may be given by: 


M{s) = 87rci(l) f z‘^Ko{b-z)dz 

Jo 


8 7rci(l) 3 
18 


I ,3 3, , 5 5, b'^s^ 

2F3 


2 2" 


2 2^ 
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- 77 I-P 2 ( { 77 }, {1, 77 }, 




log 




9 'i’{0,1+ j) 

9 2^ 


j=0 


(3 + 2 j)r 2 (i + j) 


(27) 


in which b stands for 6(1) and r(a:) and 'l>{n,x) denote the gamma function and its 
(n+1)*^ derivative (i.e, the digamma - or psi - function in the case n = 0, and polygamma 
function in the general case). Both expressions, Eq. ( p^ ) and Eq. (27) have similar perfor¬ 
mances in Mathematica. The total mass, M{oo) = ^exp6, is obtained to within 10“^, 


for s = 10. 


For the gravitational potential, using Eqs. ( 0 ) and (]^), we find the following ex¬ 
pression: 


lE'(s) = 7rci(m) s * 


,, 1 m-1, ,,, , ,, , 2m-|-l m-l-1,, , . n 

^ {}}’ {{/3s(m)}, {-,-C 2 (w) s 


2m 2m 


_ / \ ui+i ^ 2m,2 I f \ z 

= 7 rci(m)s ™ G2^2m+2 C2{m) s 


2m 2m 

{-(^ )'(^)}d } ^ 

{/3s(m)}.{-(2^).-(l5^)}, 


(28) 


Notice that, as before with the equations for v{s) and for the same reasons, in the 
case m = 4, Eqs. ( p^ and ( p^ ) above will differ from those given before in Section 4.1 
(Eqs. (|^ and (|^)). 

The gravitational potential energy is defined by: 

Vt{s) = ^ / 4>(s)dM(s) = ^ / s^i/(s)4>(s)ds (29) 

2 7o 2 7 q 

Unfortunately there seems to be no formal solution for this integral in terms of Meijer 
functions. However, by making use of a classical integral of a product of Meijer functions 


given in http://functions.wolfram.com/07.09.16.0025 one finds the following expression 
for the total potential energy: 


7rcf(m) 
“(00) = —5+37^- 


4 Cn 


■(m) 


«I ■ I (*<"■)).(-! - ->/2 - ^ 1 


lM 


' Lt. 


4 c 


2 ' (rn) 


'2m,2m+2 

2m+2,2m+2 


1- ^ - i- ,{} 


(30) 


where, besides the 2m—array of coefficients /3s(m) defined by Eq. (pSf), we have also 
defined the 2m—array: a£;(m) = —(m -I- l)/m — j3s{m). 

The velocity dispersion of a spherical system in hydrostatic equilibrium is given by: 


\s) = 


1 


M{x)v{x) 


v{s) 


dx 


(31) 





















A. Mazure & H.V. Capelato: Exact solutions... 


9 


As for the gravitational potential energy, for systems endowed with a Sersic density 
profile, this integral may be expressible in terms of Meijer functions for the case s = 0. 
We proceed similarly as to Eq. (^) to get: 

^ ({‘-i- 0 }'■ 0 " 

_ 1 TT c1{m) „2m,2m+l 

" HO) 

with aa{ra) = {2m — \)/2m — j3s{m). 

5. Conclusions 

1. We obtain analytical solutions for the de-projected de Vaucouleurs and Sersic laws 
using formal integration with Mathematica as well as for other related quantities like 
the mass or the potential, total potential energy and the central velocity dispersion. 

2. Comparisons with existing numerical estimates show very few differences, however 
analytical expressions are always much more very convenient to deal with in many 
cases. 
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Appendix A: 


The projected luminosity profile is defined by: 

fR pX — RfRe 

L{R) = 2tt 


0 


I{R')R' dR' = 2 ttIoRI 


0 


i{x')x' dx' 


(A.l) 
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l{x) = 


In terms of the dimensionless quantities, x = R/Re and l(x) = L/IqR^, this gives: 

2m7re'’ 2m7re'’ r i 

i (2m) — i (2m, bx’^) 


■ 7(2m, hx’^) = 


(A.2) 


^2m ' '' ’ ' ^2rt 

where ^{a,x) is the incomplete gamma function and F its complement (cf. Gradshteyn 
& Ryzhik, 1980). For integer values of m it can also be expressed as: 

2rmre^ 


l{x) = 


b2r. 


(2m - 1)! - exp {-hx^/^) ( (2m - 1)! + ^ (2m _ j)! 


(A.3) 


(A.4) 


From this one may find h{m) as the solutions of the equation L{Re) = Ltotj^-, where 
Ltot is the total luminosity integrated to infinity. We have (see also Ciotti 1991): 

7 ( 2 m, b) = r(2m, b) = r(2m)/2 

This is solved instantaneously using the following Mathematica commands: 

mm = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15} 
b[m_] := FindRoot[Gamma[2m, b] == Gamma[2m]/2, {b, 2m - 1/3}, 
WorkingPrecision -> 60, AccuracyGoal -> 30] 
blist = Map[b, mm]; bb = N[b /. blist, 17] 


The derived values of 6(m) are given in Table A.l , where we compare the Mathematica 
results with those from the asymptotic expansions by Ciotti & Bertin (1999). For m = 4, 
we find of course &(4) = 7.66925. 


Table A.l. The values of 6(m) 


m 

b{m) 

bcB99 

b(m) - bcB99 

1 

1.67834699001666 

1.67838865492157 

-4.16649049 10"® 

2 

3.67206074885089 

3.67206544591768 

-4.6970667 10“® 

3 

5.67016118871207 

5.67016250849902 

-1.3197869 10“® 

4 

7.6692494425008 

7.66924998466950 

-5.421687 10“’^ 

5 

9.66871461471413 

9.66871488808778 

-2.733736 10“’^ 

6 

11.6683631530448 

11.6683633097115 

-1.566667 10“’' 

7 

13.6681145993449 

13.6681146973462 

-9.80013 10"® 

8 

15.6679295443172 

15.6679296096535 

-6.53363 10"® 

9 

17.6677864177885 

17.6677864635090 

-4.57206 10"® 

10 

19.6676724233057 

19.6676724565414 

-3.32357 10"® 

11 

21.6675794898319 

21.6675795147457 

-2.49138 10"® 

12 

23.6675022752263 

23.6675022943807 

-1.91544 10"® 

13 

25.6674371029624 

25.6674371180047 

-1.50423 10"® 

14 

27.6673813599995 

27.6673813720274 

-1.20280 10"® 

15 

29.6673331382212 

29.6673331479896 

-9.7684 10"'’ 
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The luminosity so writes: 

• for m = 2, with b = b{2): 

^ [6 — {& + & by/x + ib"^ x + b^ )] (A.5) 

• for m = 4, with b = fo(4): 

|?e^[5040 - 

^-bxi (5040 + 5040 6x3 +2520 6^7^+ 840 63x3 + 

21064x +426^x3+766x5+6^x3)] (A.6) 

and so on ... 



